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Diffraction refers to the interference phenomena of propagating wave when encountering scatters, 
and has been a powerful tool for revealing the geometry. Inelastically scattered photons also carry 
phase information from internal degrees of freedom such as spin of the scatterp"*' a new ingredient 
for diffraction phenomena. Here we exploit the Fraunhofer diffraction of Stokes photons for detecting 
spin entanglement and for probing spin and motional dynamics in cold atomic ensembles. For many- 
body states where the spin-excitation number has small or maximum uncertainty, entanglement is 
witnessed by a qualitative diffraction feature, namely a sharp peak or dip in a chosen direction. 
For general states, the spin pair-correlation sum can be read out from the relative strength of 
the peak or dip over its background which, together with the expectation value and uncertainty 
of spin-excitation number, detect entanglement through the spin squeezing inequalities.^^ Spin 
precessions in Zeeman field gradient lead to displacement of the diffraction peak or dip while atomic 
motions lead to decay of its strength. These can serve as principles for vector gradiometer of fields 
and for non-demolition measurement of atomic temperature and collisional dynamics. Using a spin- 
coherent-state of TV unentangled atoms as the probe, the information of interest are encoded by 
a collectively-enhanced diffraction peak of strength ~ TV 2 , making possible gradiometer with 1/7V 
sensitivity and thermometer by a single shot. 



For atoms in a A configuration, two hyperfine ground 
states can be coupled to a common excited state by 
optical fields of different frequencies. With one of the 
transitions coupled by a laser, a spin-excitation in the 
ground state manifold can be converted into a Stokes 
photon through a spontaneous Raman transition with 
phase information preserved. This process has been ex- 
tensively explored as an interface between flying single 
photons and quantum memory based on collective spin- 
excitations in atomic ensemblesP^With a detuned laser, 
the conversion of collective spin-excitations to Stokes 
photons has the same description to the collective spon- 
taneous emission by excited two- level atomsP^The latter 
has been another topic of fundamental interest especially 
in connection with the phenomenon of super-radiance 
from very dense ensembles where atoms are close together 
compared to the optical wavelengthp2303 Recent stud- 
ies have extended this phenomenon to a new paradigm, 
namely the single photon super-radiance by ato mic e n- 
sembles prepared with single collective excitation! 14 ' 18 * ^ 

Here we exploit a sharp feature in the Fraunhofer 
diffraction of Stokes photons for detecting spin entangle- 
ment and for probing spin and motional dynamics in cold 
atomic ensembles, where interatomic distance is typically 
comparable to or larger than the optical wavelength. We 
show the sum of spin pair-correlations in a general many- 
body state can be read out from the relative strength of a 
sharp diffraction peak or dip over its background (Fig. 1). 
This pair-correlation sum together with the expectation 
value and uncertainty of spin-excitation number are suf- 
ficient for detecting e ntang lement through available spin 
squeezing inequalities ! 12 ' 13 ' expressed as pair-correlation 
sum rules. In the limit of small uncertainty of spin- 
excitation number entanglement is verified simply by see- 



ing either a diffraction peak or dip, while in the opposite 
limit of maximum uncertainty entanglement is verified 
by seeing a dip. Conserved by atomic collisions, the pair- 
correlation sum only decays in a slow timescale equal to 
the single spin dephasing time, hence this entanglement 
detection is robust and applicable for cold atom gas as 
well as atoms in optical lattices. In ensembles with many 
spin-excitations, the large photon counts from each indi- 
vidual measurement cycle can greatly facilitate the de- 
tection. Entanglement between excited two-level atoms 
can be detected in the same way from the angular distri- 
bution of spontaneously emitted photons. 

Inhomogeneous spin precessions in a Zeeman field gra- 
dient lead to displacement of the diffraction peak (dip) 
preserving its shape and relative strength. This makes 
possible calibration of inhomogencity in optical lattices 
with a precision limited only by the single spin homoge- 
neous broadening. It is also a principle for ultra-sensitive 
vector gradiometer of fields. The gradient is simultane- 
ously measured by the ~ ./V 2 pair-correlations of a TV- 
body wavefunction. Using a spin-coherent-state as the 
probe state, the gradient is encoded as the displacement 
of a diffraction peak of strength ~ TV 2 , an analog of 
the multi-slit intcrfcrometry. As a result, the gradiome- 
ter sensitivity has the 1/N scaling, an example of col- 
lectively enhanced metrology without entanglement JHHIH 
In trapped atomic gases, motional dynamics results in 
temporal decay of the peak which can be used for non- 
demolition probe of temperature and collisional interac- 
tions. 

Consider an optically thin cold atomic ensemble driven 
by a detuned laser with wavevector ko = feoz, and 
atomic motion can be considered as frozen in the du- 
ration of photon emission. Emission of a Stokes pho- 
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FIG. 1: Fraunhofer diffraction of Stokes photons from permutation-symmetric dilute ensembles. Atoms are driven 
by a detuned laser with wavevector kn = fcnz. The pair-correlation sum P in the many-atom state of interest manifests as a 
sharp diffraction peak (for positive P) or dip (for negative P) along the forward direction, with a width inversely proportional 
to the ensemble size. The value of P can be read out from the relative strength of the peak or dip to its background, provided 
that photon emission is either much faster (for cold atoms) or much slower (for hot atomic vapor) than atomic motions. This 
diffraction feature can be measured in a chosen direction with finite 6 to avoid laser photons when a Zeeman field gradient 
imprints a phase gradient on the initially permutation-symmetric ensemble (cf. Fig .3). 



ton into mode k = (k,8,ip) where k = kg is ac- 
companied by the annihilation of a spin excitation by 
J-(Ak) = Vj" , <■ lAkr t, . Ak = k k . The an- 
gular distribution of the photon emission rate is given 
by I(9,ip,t) = I s (8)I c (8,ip,t). I s is the single atom 
dipole emission pattern, a slow varying function of 8. 
I c (6,tp,t) = Tr [ J+ ( Ak) J' ( Ak) p(t) ] is the collective fac- 
tor where p(t) is the atomic density matrix. At the initial 
time of photon emission, 

/ C (0,V,O) = (N s ) + Y,e-' lA ^- rjl) (tpi) (1) 



-iAk-r, \ 12 



- ^-—l +P N 2 -N ■ 

where N s = + l)/2 is the spin-excitation num- 

ber operator. Here and hereafter (• • • ) denotes the ex- 
pectation value over p(0), the many-body state of inter- 
est. P = (J2j^j> IS thus the sum of spin pair- 
correlations. The last equal sign in Eq. holds when 
p(0) is invariant under permutation of atoms, which is 
the typical situation for atom gases. For general states 
in optical lattices, by sudden release of atoms into a 
spin-independent trappy the density matrix averaged 
over many ensemble copies will become permutation- 
symmetric after atoms lose memory of their initial po- 



sitions, with P preserved. |(X^ 7 ' e 



iAk-r, 



is a sharp 



feature which equals to N 2 along the forward direc- 
tion (8 — 0), and vastly drops to zero for 8 > 8^ = 

rnin{yS^y, ^aA^ where A and H are respectively the 
transverse and longitudinal size of the ensemble (Fig. 1). 



Thus, positive (negative) pair-correlation sum manifests 

as a sharp diffraction peak (dip), and P can be read 

out from the ratio of the peak (dip) to the background: 
i(e=o)~i(e b ) _ P 

I (8b) (N s )-P/N' 

With the spin-excitation number conserved in many 
physical processes, its expectation value (N s ) and un- 



certainty AN S = J (Nf) — (N s ) 2 are usually known a 
priori. The peak (dip) to background ratio is then suffi- 
cient to detect entanglement via available spin squeez- 
ing inequalities derived for first and second moments 
of total spinP^HUl The longitudinal component of to- 
tal spin is equivalent to the spin-excitation number: 
J z = N s — N/2, and the second moment of transverse 
components is equivalent to the pair-correlation sum: 
(J 2 ) + (J 2 ) = P + N/2. Many spin squeezing inequalities 
can thus be formulated as pair-correlation sum rules. For 
example, the optimal spin squeezing inequalities discov- 
ered in Ref. 12( become: 



(N- 



P < (N — l)AN 2 
P > -AN 2 , 

1)P > (^v s 2 ) - N(A 



(2a) 
(2b) 
(2c) 



Violation of any one of Eq. ( 2apc ) implies entanglement. 

Entanglement detection based on the above pair- 
correlation sum rules is described by the phase diagrams 
shown in Fig. 2(a-c). For eigenstates of J 2 and J z , the 
phase space characterized by the total spin quantum 
number J and the magnetic quantum number M con- 
sists of two continuous regions corresponding to diffrac- 
tion peak and dip in the forward direction, respectively. 
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FIG. 2: Entanglement detection from the diffrac- 
tion peak or dip. (a) Phase diagram in the parame- 
ter space ((N 3 ), AN a , P), being respectively the expectation 
value and uncertainty of the spin-excitation number, and the 
pair-correlation sum. States in the surrounded region are all 
entangled ones, (b) A slice of the diagram in (a) taken for 
half-spin-excitation states, i.e. (N s ) — N/2. The red, blue 
and green regions are entangled states violating inequalities 
Eq. (2a I, (2b I and pq) respectively. The grey surfaces in 



(a) and the black curves in (b) are boundaries between phys- 
ical and unphysical regions. Positive and negative sections 
of P axis have different linear scale, (c) Strength of the 
diffraction peak or dip for eigenstates of total spin J 2 and 
J z . This phase space consists of a peak region and a dip re- 
gion, violating Eq. (2a I and (2b l respectively. States violating 
Eq. pel are to the left of the dashed curve, (d) Upper (lower): 
peak (dip) to background ratio as a function of the collection 
interval r c for a half-spin-excitation state with P = 2.5iV 
(P — —0.35N), shown as the black curve. The calculation 
is for N = 4000 atoms of a 2D Gaussian distribution with 
FWHM A — 100 pm. Peak or dip (background) strength is 
evaluated at 8 — (6 = j 2 ^), shown by the blue (red) solid 
curve. Dashed curves are calculations with the multiple-light 
scattering and dipole-dipole interaction neglected. 



States in the peak (dip) region are all entangled ones 
(Eq. ([2b])). States violating Eq. (l2c 



violating Eq. (2a 



form a subset of the dip region. For a general many-body 
state, the diffraction pattern is simply the weighted av- 
erage of these patterns associated with definite J and 
M, since J2j^j> ^j^f does not mix states of different 
J and M. Qualitative criteria become possible for en- 
tanglement witness in two limits. With vanishing AN S 
seeing either a sharp diffraction peak or dip verifies entan- 
glement, while with maximum AN S seeing a dip verifies 
entanglement (Fig. 2(a-b)). On the other hand, a peak 
(dip) strength exceeding some threshold value always im- 
plies entanglement. Taking half-spin-excitation states for 
example, observing a dip to background ratio \r 



a peak to background ratio r > 



N(N-l) 
N+l 



> ± or 



ment for any possible AN S . Moreover, P and AN S can 
also quantify the entanglement depth in the vicinity of 
Dicke states!-^ 

Through these pair-correlation sum rules, the ini- 
tial diffraction pattern detects entanglement in non- 
condensed atomic ensemble regardless of shape and den- 
sity. Hereafter, we focus on dilute ensembles where inter- 
atomic distance is comparable to or larger than optical 
wavelength. Remarkably, under this condition, one can 
collect all Stokes photons, not only those initial ones, for 
measuring the pair-correlation sum in p(0). The diffrac- 
tion pattern at an arbitrary time is determined by the 
instantaneous atomic density matrix p(t) which defers 
from p(0) with spin-excitations annihilated. Neverthe- 
less, in dilute ensembles, the evolution of p under the 
spontaneous Stokes scattering is dominated by indepen- 
dent decay of spins which preserves the initial diffrac- 
tion pattern (see methods section). Corrections from 
the multiple-light scattering and dipole-dipole interac- 
tion, which need nonperturbative treatment in the dense 
limit 11SUS1 are well accounted here by perturbation theory. 
We find these weak processes only result in slow varying 
modulation of the diffraction pattern, barely changing 
the ratio of the sharp peak or dip to its neighboring back- 
ground. Namely, 



_ n(9 = 0,t c ) - n(9 b ,T c 
n(9 b ,T c ) 



P 

(N a ) - P/N 



(3) 



verifies entangle- 



where n(9, ip, t c ) = 5Q J Q C dtl(9, (p, t) is the number of 
photons emitted into an infinitesimal solid angle in the 
collection interval < t < r c . For a dilute ensemble 
Eq. ([3| holds for arbitrarily large r c (cf. Fig. 2(d)). This 
is greatly helpful in ensembles with large (N s ), since large 
number of photons can be collected in each measurement 
cycle to facilitate the entanglement detection. 

We estimate the range of applicability of our treat- 
ment. The dilute condition is satisfied by typical cold 
atom gases of a density 10 10 — 10 12 cm -3 or by atoms 
in optical lattice. The duration of Stokes photon emis- 
sion is of a timescale T _1 = (§^)~ 2 rQ 1 , where CIl is the 
laser Rabi frequency and A the detuning. The excited 
state decay rate To > 30 ps^ 1 for typical alkali atomsP^ 
Taking (§j)~ 2 ~ 40, all Stokes photons are emitted in a 
timescale < /is. For cold atom gases with a temper- 
ature of 1 — 100 /iK, the average velocity is 0.01 — 0.1 m/s. 
Hence atoms can only travel 10 — 100 nm in the duration 
of r _1 which is negligible as compared to the light wave- 
length. Moreover, we find that pair-correlation sum of a 
dilute hot atomic vapor can be measured in the same way 
if Stokes photon emission is controlled to be much slower 
than atomic motions (see supplementary material). 

Under free evolution, the pair-correlation changes as 
Tt[a+a7p( T )] = e ^*o-^T)-2 7 T Tr ^+ (5 .- p ^ where % 

is the Zeeman frequency and 7 the homogeneous dephas- 
ing rate of an individual spin. The pair-correlation sum 
thus decays only at the single spin dephasing rate. There- 
fore entanglement in p(0) can be reliably detected from 
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FIG. 3: Diffraction patterns (lower parts) from 
atomic ensemble in different spin configurations (up- 
per), (a) Spin-coherent-state with in-plane polarization, (b) 
Evolution in a Zeeman field gradient imprints a phase gra- 
dient of spins, resulting in a displacement of the diffraction 
peak which can be a principle of gradiometer. (c) Atomic mo- 
tions diminish the spin polarization, resulting in decay of the 
displaced peak. This can be a principle for non-demolition 
measurement of atomic temperature and collisional interac- 
tions. 



f ; 

^ peak decay: (3,0) 2 (Av^ 



the dephased state p(r) as long as r <C 7 1 , even when 
the fidelity is exponentially small with iVP^ 

Spatial inhomogeneity of external fields leads to posi- 
tion dependent Zeeman frequency r)(r) and hence inho- 
mogeneous precessions of spins. If the size of the ensem- 
ble is small compared to the variation length scale of the 
field, the dominating term is the gradient: r/(r) = r • Vrj. 
For an ensemble initially in a permutation-symmetric 
state, after an interval tq with frozen motion in the 
Zeeman field gradient, the diffraction pattern becomes 

(E .e- i ( Ak - ToV ")-^}| 2 . Wefo- 



I C — (NgJ J y_ 1 T N 2_ N I 

cus on situations where d z rj is either zero or not picked 
up by atomic ensembles of a quasi-2D geometry in x — y 
plane. The in-plane gradient simply results in a displace- 
ment of the sharp diffraction peak or dip, preserving its 
strength and shape (Fig. 3). This has several significant 
consequences. First, by evolution in an external field of 
known gradient, entanglement can be detected by mea- 
suring the peak or dip along a chosen direction with fi- 
nite 9, such that detectors do not pick up laser photons. 
Second, the displacement measures the vector value of 
unknown gradient. This can be used to calibrate inho- 
mogeneous broadening in optical lattices. The Zeeman 
field gradient can come from the gradient of magnetic 
field, static electric field via dc Stark effect, and light 
field via ac Stark effect!^ The physics discussed here can 
thus be used as a principle of vector gradiometer of these 
fields with a spatial resolution equal to the ensemble size 
A. 

An ideal probe state is simply the spin-coherent-state 
of N unentangled atoms with in-plane polarization. The 
gradient is then probed simultaneously by the ~ N 2 clas- 
sical pair-correlations, and its vector value is encoded 
as the displacement of a diffraction peak with strength 
~ A 2 . This A^ 2 scaling of the peak strength suggests an 
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FIG. 4: Zeeman field gradiometer using a ID atomic 
ensemble, (a) and (b) Schematic of the setup. The en- 
semble and a ID CCD array are placed respectively on the 
two common focus lines of a group of elliptical cylinder mir- 
rors. The smallest ensemble can be an atom-pair, giving a 
direct analog of the double-slit interferometry. (c) Signal 
of N/2 probes using atom-pairs (upper), and single probe 
using TV atoms collectively (lower), with each atom placed 
randomly on the focus line according to a Gaussian distri- 
bution with FWHM A. The peak strength of the lower is 
enhanced by a factor of N. (d) Gradiometer sensitivity at a 
spatial resolution A = 1 mm with a resource of N unentan- 
gled atoms. The diffraction based gradiometer using N atoms 
collectively (N/2 atom-pairs independently) has a sensitivity 
of the 1/N (1/yN) scaling, shown by the solid (dashed) black 
line. Sensitivity of flying atom Mach-Zehnder interferometry 
(MZI) gradiometer is shown for referenced The probe time 
tq — -A- = 0.1 s for the blue line, limited by a finite velocity 
uncertainty An = 1 cm/s, while to = 1 s for all other lines, 
limited only by the single spin dephasing time. 



enhancement of the sensitivity by using large group of 
atoms collectively, as compared to uncorrelated probes 
using small groups. In a ID geometry (Fig. 4), a single 
collective probe using A~ atoms can achieve a sensitivity 
of ~ ir — 7=f, which goes beyond the standard quan- 

' for N/2 independent probes 



/k A 

turn limit of ~ 



•/N ToVkoA 

using atom-pairs. This collectively-enhanced sensitivity 
with 1/N scaling is valid in the dilute regime A^ < k A. 
Beyond this regime, the multiple light scattering can not 
be treated perturbatively and its effect will eventually 
renormalize the peak strength to the A^ scaling. By trap- 
ping atoms in optical lattices, the probe time To can be as 
long as the single spin homogeneous dephasing time, in 
the order of second or longerPThis diffraction based gra- 
diometer using stationary atoms is immune to collective 
noises and uncertainty in atomic positions. In compari- 
son, the gradiometer based on Mach-Zehnder interferom- 
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etry of flying atoms in atomic fountain h as th e standard 
quantum limited sensitivity of ~ h= | 2Q | 30 | ^ut ^ e m _ 
evitable velocity uncertainty sets a tighter upper bound 
for To dependant on the spatial resolution A (see section 
VD in supplementary material). 

Finally, we present a principle for non-demolition 
probe of atomic motions and temperature in trapped cold 
atom gases (see Fig. 3). A phase gradient is first im- 
printed on the spin-coherent-state!^ After a waiting time 
Ti , we turn on the laser and measure the diffraction pat- 
tern of Stokes photons. Atomic motions in the interval 
Ti will diminish the spin polarization, resulting in decay 
of the displaced diffraction peakP For |V(/>| 2 (Ar 2 ) <C f, 
the peak strength is i^ e -|V0| 2 <Ar 2 )/3 ; ^r 2 ) being the 
mean square displacement of atoms. For short t\ when 
Ar is small compared to the interatomic distance d, 
(Ar 2 ) = 2^pT 2 . Thus, by preparing a large phase gra- 
dient |V0| ~ \, this short time motion of atoms can be 
probed and the atomic temperature can be read out from 
the decay of the peak. Smaller |V(/>| allows the probe 
of long time motion which will eventually crossover to 
the diffusive regime by atom collisions. t\ is upper lim- 
ited by the spin dephasing time, which is long enough 
for observing the entire crossover behavior from ballis- 
tic to diffusive motions, providing information about the 
collisional interactions in trapped cold atom gases. The 
collectively enhanced peak strength of ~ N 2 provides 
sufficient signal-to-noise ratio for determining (Ar 2 ) at a 
given Ti by a single shot measurement. 

Methods 

Light-atom interaction. Consider atoms of a A level 
configuration where two atomic ground states \g) and \s) 
can be optically coupled to a common excited state |e) 
(Fig. 1 inset). With a coherent laser of Rabi frequency 
JIl, detuning A, and wavevector ko = koz coupling the 
\s) to |e) transition, an atom can go from state \s) to |<?) 
by emitting a Stokes photon into the vacuum. When A is 
much larger than f2 l and the excited state homogeneous 
line width r , |e) can be adiabatically eliminated, leading 
to the effective light-atom coupling in the electric-dipole 
and rotating wave approximation: 

H = ^huikalau + ^Ezaj 

k j 



where g^ — ^-y^T^ek • £t, t?k and /x being respec- 
tively the unit polarization vector and the single atom 
dipole. The ac Stark shift has been grouped into E x , 
and the pseudospin operators are defined in the ground 



state manifold: cr- = |<7)j(s| and o* 



\9)M 



We assume anti-Stokes scattering is cither forbidden by 
the polarization selection rule or suppressed by the much 
larger detuning when E z 3> A. 

The same Hamiltonian can also describe spontaneous 
emission of excited two-level atoms where the excited 



state |e) decays to ground state \g) by emitting a pho- 
ton into the vacuum. In such case, g^ = 



e k ■ At 



and the pseudospin operators cr- = e j \g)j(e\ and 
ct| = |e)j(e| — \g)j(g\ are defined for the j-th atom lo- 
cated at Yj. 

Perturbative solution of master equation. Under 
the above light-atom interaction, the evolution of the 
atomic density matrix is described by the Lindblad mas- 
ter equation in the Born-Markov approximation: 



dp 
dt 

C oP (t) 
C lP (t) 



C oP (t) + C lP (t) 

v r 



2 



,sin(fco|rj-r 3 ,|) 



where T jf = r °";!7 ] ^"' and G iV = T 



s(fc |rj — r ■/ 1) 



03 fc |r.;-iy| 

describe respectively the mul tiple light scattering and 
dipole-dipole interaction! 14 ! 17 ! In the study of super- 
radiance phenomena in very dense atomic ensembles, 
these effects must be accounted non-perturbativelyP^ 
Tjji and G$y drop fast with distance. In dilute atomic 
ensembles where the atom-atom distance is comparable 
or larger than the photon wavelength, the atomic evo- 
lution can be solved perturbatively. Using the Laplace 
transform w{z) — dte~ zt p(t), we have 



w(z) 



1 



-p(0). 



For ko\rj — Vji\ > 2ir, L\ is small compared to Cq, and 
we make a perturbative expansion z _ c ^_ Cl = z -c + 
z 2c C\ + • • ■ ■ By inverse Laplace transform we can 
get the solution of the atomic density matrix p^ n ' keeping 
up to the n-th order effects of C\. 

For the zeroth order solution p^'{t) — e Cot p(0), we 
find Tr[J+(k)J-(k)/j( )(i)] = e- r *Tr[J+(k)J-(k)p(0)], 
i.e. the initial diffraction pattern is preserved for all time. 
Comparisons with exact solution of master equation for a 
chain of 12 atoms show that the perturbation expansion 
converges fast for k \rj — Vy \ > 2ir, and the effects of the 
multiple light scattering and dipole-dipole interaction are 
well accounted by keeping only the first order effects of 
the L\ term: 



p(V(t)=p(.o)(t)+ / dTe CoT \c lP {0) {t-T) 



The L\ term leads to slow varying modulation of the 
diffraction pattern, which barely changes the ratio of the 
sharp peak (dip) to its neighboring background. Details 
on this modulation and the convergence check for the 
perturbative solutions can be found in the supplementary 
material. 
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